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Abstract. The recent detection of high energy 7-rays 
coming from supernova remnants and active galactic nu- 
clei has revived interest in the diffusive shock acceleration 
of electrons. In the present paper we examine the basis 
of the so-called "box" model for particle acceleration and 
present a more physical version of it. Using this we deter- 
mine simple criteria for the conditions under which "pile- 
ups" can occur in shock accelerated electron spectra sub- 
ject to synchrotron or inverse Compton losses (the latter 
in the Thompson limit) . An extension to include nonlinear 
effects is proposed. 
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The EGRET detection aboard the Compton Gamma Ray 
Obse rvatory of at least two supernova remnants (Esposito 
|1996|) and m ore than fifty active galactic nuclei (Thomson 



et al. 1995 ) has given strong evidence of particle acceler- 



ation in these objects. This evidence is strengthened even 
more by the detection of SN 1006 (Tanimori 1998 ) and the 
BL Lac objec ts Mk n 421 (Punch et al. |1992D and Mkn 501 
(Quinn et al. 1996) by ground based Cherenkov detectors 
at TeV energies. 

A particularly attractive mechanism for producing the 
required radiating high energy particles is the diffusive 
shock acceleration scheme, which has already been put for- 
ward to predict TeV radiation from supernova remnants 
(Drury, Aharonian & Volk |1994 Mastichiadis |1996D or 



explain the observed flaring behaviour in X-rays and TeV 
7-rays fr om a ctive galactic nuclei (Kirk, Rieger & Mas- 
tichiadis 1998). This scheme was originally proposed as the 
mechanism responsible for producing the nuclear cosmic 
ray component in sh ock w aves asso ciated with supernova 
remnants (Krymsky |l977 , Axford 1981). Based on this 
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picture many authors (Bogdan & Volk 1983, Moraal 
Axford 



1983 



Lagage & Cesarsky 1983 



Volk & Biermann 



1988 



Schlickeiser 



1992 



1984, 



Protheroe & 



BaU & Kirk 

Stanev 1998 ) have used, under various guises, a simplified 
but physically intuitive treatment of shock acceleration, 
sometimes referred to as a "box" model. 

In this paper we examine the underlying assumptions 
of the "box" model (§2) and we present an alternative 
more physical version of it (§3). We then include syn- 
chrotron and inverse Compton losses as a means of spectal 
modification and we determine the conditions under which 
"pile-ups" can occur in shock accelerated spectra (§4). The 
"box" model can also be extended to include the nonlin- 
ear effect of the particle pressure on the background flow 
(§5). 

2. The "box" model of diffusive shock 
acceleration 

The main features of the "box" model, as presented in 
the literature (see references above) and exemplified by 
Protheroe and Stanev (1998)) can be summarised as fol- 
lows. The particles being accelerated (and thus "inside 
the box") have differential energy spectrum N{E) and 
are gaining energy at rate r^ccE but simultaneously es- 
cape from the acceleration box at rate Tcsc- Conservation 
of particles then requires 



r.N 



(1) 



where Q{E) is a source term combining advection of par- 
ticles into the box and direct injection inside the box. 

In essence this approach tries to reduce the entire ac- 
celeration physics to a "black box" characterised simply 
by just two rates, fesc and race- These rates have of course 
to be taken from more detailed theories of shock accelera- 



tion (eg Drury 1991). A minor reformulation of the above 
equation into characteristic form. 
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is useful in revealing the character of the description. This 
is equivalent to the ordinary differential equation, 



dN 
~dt 



N 



E 



dE 



on the family of characteristic curves described by 
at 

giving the formal solution, 

N{E,t)^ f Q {t',E') 
Jo 

Br, 



(3) 



(4) 



If we measure p, the magnitude of the particle's momen- 
tum, in this frame, it is not changed by the scattering 
and the angular distribution is driven to being very close 
to isotropic. However if a particle crosses a shock front, 
where the bulk plasma velocity changes abruptly, then 
the reference frame used to measure p changes and thus p 
itself changes slightly. If we have an almost isotropic dis- 
tribution f(p) at the shock front where the frame velocity 
changes from Ui to U2, then it is easy to calculate that 
there is a flux of particles upwards in momentum associ- 
ated with the shock crossings of 
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The number of particles at energy E and time t in the 
"box" is given simply by an exponentially weighted in- 
tegral over the injection rate at earlier times and lower 
energies. Of particular interest is the steady solution at 
energies above those where injection is occuring which is 
easily seen to be a power-law with exponent 
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(6) 



At first sight (to one familiar with shock acceleration 
theory) it appears odd that the exponent depends not 
just on the ratio of Tesc to race but also on the energy de- 
pendence of race- However, as remarked by Protheroe and 
Stanev, the physically important quantity is not the spec- 
trum of particles inside the fictitious acceleration "box" 
but the escaping flux of accelerated particles VcscN and 
this is a power-law of exponent 



dlnjrescN) 
d\nE 
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Thus provided the ratio of race to rose is fixed, the power- 
law exponent of the spectrum of accelerated particles es- 
caping from the accelerator is determined only by this 
ratio whatever the energy dependence of the two rates. 

3. Physical interpretation of the box model 

We prefer a very similar, but more physical, picture of 
shock acceleration which has the advantage of being more 
closely linked to the conventional theory. For this reason 
we also choose to work in terms of particle momentum 
p and the distribution function f{p) rather than E and 
N{E). 

The fundamental assumption of diffusive shock accel- 
eration theory is that the charged particles being accel- 
erated are scattered by magnetic structures advected by 
the bulk plasma flow and that, at least to a first approxi- 
mation, in a frame moving with these structures the scat- 
tering changes the direction of a particle's motion, but 
not the magnitude of its velocity, energy or momentum. 



where n is the unit shock normal and the integration is 
over all directions of the velocity vector v. Notice that this 
flux is localised in space at the shock front and is strictly 
positive for a compressive shock structure. 

This spatially localised flux in momentum space is the 
essential mechanism of shock acceleration and in our de- 
scription replaces the acceleration rate race- The other key 
element of course is the loss of particles from the shock by 
advection downstream. We note that the particles inter- 
acting with the shock are those located within about one 
diffusion length of the shock. Particles penetrate upstream 
a distance of order Li = n • Ki • n/n • Ui where K is the 
diffusion tensor and the probability of a downstream par- 
ticle returning to the shock decreases exponentially with a 
scale length of L2 = n • K2 • n/n • U2. Thus in our picture 
we have an energy dependent acceleration region extend- 
ing a distance Li upstream and L2 downstream. The total 
size of the box is then L{p) = Li{p) + L2{p). Particles are 
swept out of this region by the downstream flow at a bulk 
velocity n • U2 . 

Conservation of particles then leads to the following 
approximate description of the acceleration, 



d (9$ 

- [4V/i] + ^ = Q - n . U24V/, 



(9) 



that is the time rate of change of the number of parti- 
cles involved in the acceleration at momentum p plus the 
divergence in the accelerated momentum flux equals the 
source minus the flux carried out of the back of the region 
by the downstream flow. The main approximation here 
is the assumption that the same f{p, t) can be used in all 
three terms where it occurs. In fact in the acceleration flux 
it is the local distribution at the shock front, in the total 
number it is a volume averaged value, and in the loss term 
it is the downstream distribution which matters. Diffusion 
theory shows that in the steady state all three are equal, 
but this need not be the case in more elaborate transport 



models (Kirk, Duffy & Gallant |1 9961 ). 

Substituting for $ and simplifying we get the equation 

i| + „.U./,p) + i„^(U.-U,)p|.j|^ ,10) 
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/Advection In 



of the "box" or region normally increases with energy this 
also gives an additional loss process because particles can 
now fall through the back of the "box" as well as being ad- 
vected out of it (see Fig. 1). Note that particles which fall 
through the front of the box are advected back into the 
acceleration region and thus this process does not work 
upstream. 

If the loss rate is p — —ap^ the basic equation becomes 



Q - C/24V/(P) - ^wap'^fip)^ 



Fig. 1. Sketch of the acceleration region or "box" indicat- 
ing the particle fluxes. The dashed lines indicate the front 
and back edges of the acceleration region. 

which is our version of the "Box" equation. Note that this, 
as is readily seen, gives the well known standard results 
for the steady-state spectrum and the acceleration time- 
scale. In fact our description is mathematically equivalent 
to that of Protheroe and Stanev as is easily seen by noting 
that 



n ■ (Ui - U2) 
3L 



N = Airp^fL-iU) 



However our version has more physical content, in partic- 
ular the two rates are derived and not inserted by hand. 
It is also important to note that in our picture the size of 
the "box" depends on the particle energy. 

4. Inclusion of additional loss processes 

In itself the "box" model would be of little interest be- 
yond providing a simple "derivation" of the acceleration 
time scale. Its main interest is as a potential tool for in- 
vestigating the effect of additional loss processes on shock 
acceleration spectra. One of the first such studies was that 



of Webb, Drury and Biermann (1984) where the impor- 
tant question of the effect of synchrotron losses was in- 



vestigated (see also Bregman et al. 1981). An interesting 
question is whether or not a "pile-up" occurs in the ac- 
celerated particle spectrum at the energy where the syn- 
chrotron losses balance the acceleration. Webb, Drury and 
Biermann (1984) found that pile-ups only occured if the 
spectrum in the absence of synchrotron losses (or equiva- 
lently at low energies where the synchrotron losses are in- 
significant) was harder than / oc p~'^. However Protheroe 
and Stanev obtain pile-ups for spectra as soft as / cx p^^-^. 

It is relatively straightforward to include losses of the 
synchrotron or inverse Compton type (Thomson regime) 
in the model. These generate a downward flux in momen- 
tum space, but one which is distributed throughout the 
acceleration region. Combined with the fact that the size 



(12) 



This equation is easily generalised to the case of different 
loss rates upstream and downstream. Simplifying equation 
(12) gives 
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Note that for convenience we have dropped the explicit 
vector (and tensor) notation; all non-scalar quantities are 
to be interpreted as normal components, that is C/2 means 
n • U2 etc. Note also that our model differs from that of 
Protheroe and Stanev in that they do not allow for the 
extra loss process resulting from the energy dependence 
of the "box" size. 

In the steady state and away from the source region 
this gives immediately the remarkably simple result for 
the logarithmic slope of the spectrum. 



dlnf 
dlnp 



2 dLi 

Ui — AapL — ap — — 
, dp_ 

Ui — U2 — 3apL 



(14) 



Note that at small values of p we recover the standard 
result, that the power-law exponent is —3Ui/{Ui — C/2). 

Under normal circumstances both Li and L2 are 
monotonically increasing functions of p. Thus both the 
numerator and denominator of the above expression, re- 
garded as functions of p, have single zeroes at which they 
change sign. The denominator goes to zero at the critical 
momentum 



P 



U1-U2 
3aL 



(15) 



where the losses exactly balance the acceleration. If the 
numerator at this point is negative, the slope goes to —00 
and there is no pile-up. However the slope goes to +00 
and a pile-up occurs if 
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In the early analytic work of Webb et al the diffusion 
coefficient was taken to be constant, so that dLi/dp = 



4 
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log(f)+qlog(p) 




Fig. 2. Energy spectra for a momentum dependent diffu- 
sion coefficient k oc and a compression ratio of r = 4 
where q = 3r/ (r — 1) 
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Fig. 3. As in figure 2 but with r = 3. 

and this condition reduces to Ui > AU2 in agreement with 
their results. However if, as in the work of Protheroe and 
Stanev, the diffusion coefficient is an increasing function 
of energy or momentum, the condition becomes less re- 
strictive. For a power-law dependence of the form K (x 
the condition for a pile-up to occur reduces to 



Ui -AU2 + 5 {Ui - U2) 



Li + L2 



> 



(17) 



(The equivalent criterion for the model used by Protheroe 
and Stanev is slightly different, namely 



because of their neglect of the additional loss process.) 

For the case where L1/L2 = U2/U1 and with 5=1 
this condition predicts that shocks with compression ratios 
greater than about r = 3.45 will produce pile-ups while 
weaker shocks will not. fn Figures 1 and 2 we plot the 
particle spectra up to p* for a range of values of 6 and 
with r = 4 and r = 3 respectively. 

Thus there is no contradiction between the (exact) re- 
sults of Webb et al and those of Protheroe and Stanev; 
the apparent differences can be attributed to the energy 
dependence of the diffusion coefficient. Indeed, looking at 
the results presented by Protheroe and Stanev, it is clear 
that the pile-ups they obtain are less pronounced for those 
cases with a weaker energy dependence. 

5. Nonlinear effects 

At the phenomenological and simplified level of the 
"box" models it is possible to allow for nonlinear ef- 
fects by replacing the upstream velocity with an effec- 
tive momentum-dependent velocity Ui{p), reflecting the 
existence of an extended upstream shock precursor region 
sampled on different length scales by particles of differ- 
ent energies. Higher energy particles, with larger diffusion 
length scales, sample more of the shock transition and 
have larger effective values of Ui{p); thus Ui{p) must be 
a monotonically increasing function of p. Repeating the 
above analysis with a momentum-dependent Ui the loga- 
rithmic slope of the spectrum is in this case 



Ui — AapL 



p dUi 
3 dp 



ap 



dlnf _ 

dlnp " Ui — U2 — 3apL 

with a pile-up criterion of. 



) dLi 
dp 



(19) 



Ui {p) - W2 -p^ +3qi/ '^^^ 



dUi 
dp 



dp 



>0 at p = p*{20) 



Ui -4U2 + 6 (t/i -U2)>0 



(18) 



We see that whether or not the nonlinear effects assist the 

formation of pile-ups depends critically on how fast they 
make the effective upstream velocity vary as a function of 
p. By making Ui{p*) larger they make it easier for pile-ups 
to occur. On the other hand, if the variation is more rapid 
than Ui (X p, the derivative term dominates and inhibits 
the formation of pile-ups. 

In most cases the shock modification will be produced 
by the reaction of accelerated ions, and the electrons can 
be treated as test-particles with a prescribed ?7i(p). How- 
ever in a pair plasma, or if one applies the "box" model 
to the ions themselves, the effective upstream velocity has 
to be related to the pressure of the accelerated particles 
in a self-consistent way. We require in the "box" model a 
condition which describes the reaction of the accelerated 
particles on the flow. Throughout the upstream precursor 
and in the steady case both the mass flux, A = pU, and the 
momentum flux, AU + Pc are conserved. Here Pc is the 
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pressure contained in energetic particles and the gas pres- 
sure is assumed to be negligible upstream. At a distance 
Li{p) upstream only particles with momenta greater than 
p remain in the acceleration region. This suggests that in 
the "box" model the reaction of the particles on the flow 
is described by the momentum flux conservation law 



constant 



(21) 



where Pmax is the highest momentum particle in the sys- 
tem and V is the particle velocity corresponding to mo- 
mentum p. Differentiating with respect to p gives 



A^=4np'fip)Pl. 

dp f jyfi 2, 



(22) 



With no losses and for Ui{p) ^ U2 we can now recover 
Malkov's spectral universality result for strong modified 
shocks (Malkov, |1998D . In the hmit of U2 = and a = 
the conservation equation reduces to the requirement than 
the upward flux in momentum space be constant (equation 

(9)), 

<i> = ^/b)f/ib)-*o- (23) 



When combined with equation (22) this gives 



(24) 



dUi _ $0 _ ^0 dT 
dp ~ A ~ A dp 

where we have used the elementary result from relativistic 
kinematics that the particle velocity is the derivative of 
the kinetic energy T with respect to momentum. Integrat- 
ing for relativistic particles, T = pc, we get the fundamen- 
tal self-similar asymptotic solution found by Malkov, 



Ui = 



2c$, 



A 



1/2 



2c 



-P 



-3.5 



(25) 



If the electrons are test-particles in a shock strongly mod- 
ified by proton acceleration, and if the Malkov scaling 
C/i oc p^/"^ holds even approximately, then equation (20) 
predicts that a strong synchrotron pile-up appears in- 
evitable. 

It is perhaps worth remarking on some peculiarities of 
Malkov's solution. Formally it has U2 = 0, all the kinetic 
energy dissipated in the "shock" is used in generating the 
upwards flux in momentum space $ and there is no down- 
stream advection. It is not clear that a stationary solution 
exists in this case. The problem is that as C/2 — *■ so 
L2 ^ 00 \i a. diffusion model is used for the downstream 
propagation. The solution appears to require some form 
of impenetrable reflecting barrier a finite distance down- 
stream if it is to be realised in finite time. Also, although 
the accelerated particle spectrum at the shock is a uni- 
versal power law, none of these particles escape from the 
shock region. From a distance the shock appears as an al- 
most monoenergetic source at whatever maximum energy 
the particles reach before escaping from the system. 



The case of a synchrotron limited shock in a pure pair 
plasma is also interesting. Here the upper cut-off is de- 
termined not by a free escape boundary condition but by 
the synchrotron losses. If most of the energy dissipated 
in the shock is radiated this way, the shock will be very 
compressive and the downstream velocity U2 negligible 
compared to [/i. The same caveats about time scales ap- 
ply as to Malkov's solution, but again we can, at least as a 
gedanken experiment, consider a cold pair plasma hitting 
an impenetrable and immovable boundary. In this case, 
if there is a steady solution, the upward flux due to the 
acceleration must exactly balance the synchrotron losses 
at all energies. In general it appears impossible to satisfy 
both this condition and the momentum balance condition 
for p < p* unless the diffusion coefficient has an artificially 
strong momentum dependence. However a solution exists 
corresponding, in the box model, to a Dirac distribution 
at the critical momentum p*. This steady population of 
high energy electrons has enough pressure to decelerate 
the incoming plasma to zero velocity and radiates away 
all the absorbed energy as synchrotron radiation. This ex- 
treme form of pile-up may be of interest as a means of very 
efficiently converting the bulk kinetic energy of a cold pair 
plasma into soft gamma-rays. 

6. Conclusion 

A major defect of all "box" models is the basic assump- 
tion that all particles gain and loose energy at exactly the 
same rate. It is clear physically that there are very large 
fluctuations in the amount of time particles spend in the 
upstream and downstream regions between shock cross- 
ings, and thus correspondingly large fluctuations in the 
amount of energy lost. The effect of these variations will 
be to smear out the artificially sharp pile-ups predicted by 
the simple "box" models. However our results are based 
simply on the scaling with energy of the various gain and 
loss processes together with the size of the acceleration 
region. Thus they should be relatively robust and we ex- 
pect that even if there is no sharp spike, the spectrum will 
show local enhancements over what it would have been in 
the absence of the synchrotron or IC losses in those cases 
where our criterion is satisfied. 
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